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1 Introduction 



A relevant class of two-dimensional integrable models are those based on the representations of a 
given quantum affine algebra U q [G] [lj. The quantum group framework permits us to generate solu- 
tions of the Yang-Baxter equation which can be seen as the Boltzmann weights of two-dimensional 
vertex models [2]. The standard paradigm is the spin-s representation of the U q [SU(2)} algebra 
for generic values of q [3] which leads us to an iV-state extension of the six- vertex model pfl, [5]. 
The corresponding one- dimensional spin chain turns out to be an integrable generalization of the 
Heisenberg model with arbitrary spin s = (N — l)/2 |6J. 

The previous example by no means exhaust the possibility of constructing vertex models within 
the quantum U q [SU{2)] algebra. It still remains the representations for non-generic values of q 
which are known to be quite different from the one with arbitrary q [7j. In fact, a new family of N- 
state vertex model can be generated by exploring non-cyclic U q [SU(2)] representations when q is a 
root of unity. The respective .R-matrix depends on both an extra continuous variable besides the 
spectral parameter and on a discrete variable characterizing the possible roots of unity branches. 
We remark that this .R-matrix has been previously discussed in the literature in distinct contexts 
such as on the realm of new braid matrices with extra color variables [Sj, in the Baxterization of 
braids for the special cases of N = 2, 3 [9] and in the quantum group framework [TOj ITT] . 

The Bethe ansatz solution of such iV-vertex model based on roots of unity has so far been 
restricted to the case of periodic boundary conditions [T2T [TB] . Results for fixed boundary con- 
ditions are concentrated on the calculations of the partition function with certain domain wall 
boundary conditions [H]. The purpose of this paper is to start the study of this class of vertex 
models with open integrable boundary terms. We shall here consider the simplest case of diago- 
nal boundary conditions containing a free parameter. We find that the dependence of the Bethe 
equations on the continuous parameter characterizing the roots of unity representation appears 
only in the pseudomomenta function. By way of contrast the respective many-body phase-shifts 
are functions determined by the discrete roots of unity branches and the boundary parameters. 

The outline of this paper is as follows. In section 2 we review the structure of the bulk .R-matrix 
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and discuss properties that are useful to build up boundary K-matrices. In section 3 we discuss 
the most general diagonal solution of the reflection equations for both left and right boundaries. 
We also present the explicit expressions of the corresponding one-dimensional open spin chains for 
N = 2 and N = 3. In section 4 we adapt the algebraic Bethe ansatz construction of [!15J in order 
to derive the eigenvalues of the double-row transfer matrix and the respective Bethe equations. 
Our conclusions are summarized in section 5. In Appendices A and B we present technical details 
concerning the form of projectors and i?-matrix amplitudes. 



2 The i?-matrix properties 

The corresponding i?-matrix of the A-state vertex model based on the U q [SU(2)] algebra at roots 
of unity has been previously discussed by several authors [3J [91 [10l [11]. We shall follow here 
the recent presentation given by us in terms of the Baxterization approach [13] . For practical 
computations it is convenient to write the i?-matrix introducing an auxiliary operator R(X), 

R 12 (X) = P 12 R 12 (X), (1) 

where P\ 2 denotes the C N <8> C N operator. 

The matrix R(X) can be expressed by means of the following linear combination of projectors, 

^)=/KA)En^K^Mfi(7.*), (2) 
i=i j=i smh[ % + 27 - A] 

where k is an integer coprime to iV and the overall normalization p(X) is chosen, 

p(A)=n 1 sinh[-A + n+ ^~ 1) ]. (3) 

i=i 

The i?-matrix ([U13]) is characterized by the continuous variable 7 and the discrete index k 
parameterizing the non-cyclic U q [SU(2)] representation at the roots of unity q = e^. The 
projectors Pi(X,k) can be written in terms of the corresponding braid representation [8] whose 
expressions can be found in [13J. To make this paper self-consistent we have summarized the main 
formula for Pj(7, k) in Appendix A. In order to investigate integrable open boundary conditions it 



is convenient to bring the bulk i?-matrix (jT}{3]) to its most possible symmetrical form. As usual this 
is accomplished by performing a transformation in i?-matrix (jT]{3]) that preserves the Yang-Baxter 
equation. We find that the suitable spectral parameter dependent transformation is, 



R 12 (\) = V 1 (X)R 12 (X)V l ~ \X). 
where the gauge matrix V(X) is diagonal and it is given by 

N 



(4) 



(5) 



a=l 



while e a> b denotes standard N x N Weyl matrices. 

Let us now discuss the structure of the R- matrix .R 12 (A) for few values of N. The simplest case 
N = 2 turns out to be directly related to a six- vertex model satisfying the free-fermion condition, 
namely 



' sinh [27 + A] 
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sinh [A] sinh [ij] 

sinh [27] sinh [A] 




\ 



(6) 



sinh[z7 — A] j 

The i?-matrix (TTJIH]) for > 3 gives origin to novel integrable vertex models. For the simplest 
case N = 3 one has a nineteen-vertex model and the structure of the corresponding .R-matrix is, 
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The respective Boltzmann weights are given by, 

ITT k 

a±(X) = sinh[z7 ± A] sinh[z7 H — — ± A] 

o 

17T k 

b + (X) = sinh[A] sinh[«7 H — - — h A], b-(X) = sinh[A] sinh[z7 — A] 

iTT k %ixk 

c + (A) = sinh[z7] sinh[z7 H h A], c_(A) = sinh[z7 — A] sinh[z7 H ] 

3 3 

/ / %'TT J^j ITT f}j 

d(X) = e k sinh[A]y sinh[z7] sinh[z7 H — — ], e(A) = sinh[z7] sinh[z7 H — — ] 

V o o 

%7T ki %ixk Z7T ' Aj 

/(A) = sinh[A] sinh[A H — — ], g(X) = sinh[z7] sinh[z7 H — — ] + sinh[A — ] sinh[A] (8) 

o o o 

where the phase Ek = exp[z7r(A; — 1)]. 

Considering Eqs.([TJl5]) and the expressions for the projectors ( 1A.lHA.4l) one is able to compute 
the i?-matrix elements of R±2 (A) on the Weyl basis, 

JV 

Ri2(X)= Yl ^(A)e„, c ®e M (9) 

a,b,c,d=l 

for rather high values of A" with moderate computational effort. As further examples we have 
exhibited in Appendix B the non-null Boltzmann weights of the R12 (A) for N = 4 and N = 5. 

We now turn our attention to describing the symmetry properties of R12 (A). Besides satisfying 
the unitarity property it is invariant by both temporal and parity symmetry, namely 

Ri2(X)R 21 (-X) = p{X)p{-X)I N ® I N (10) 
Pi2^i 2 (A)Pi 2 = Ri 2 {X) (11) 
R 12 (Xy^ = R 12 (X) (12) 

where In is the N x N identity matrix and £3 denotes the transposition on the z'-th space. 

We also note that for arbitrary 7 the matrix elements -R^(A) of ^12 (A) are not self-conjugated 
under the charge invariant*]- As a consequence of that the .R 12 (A) does not satisfy standard 
crossing properties but only an analog of that which is given by, 

i?i 2 (A)' 2 i? 21 (-A + ^)' 2 = Pl (X)I N ® I N (13) 



1 Recall that the charge symmetry relates the weights R% b (\) and R N +\_a%^\_ b (X) 



where pi(A) = Yl fji sinh 2 [A + v ,. 

In next section we shall use the properties (11011131) to obtain the right .KT-matrix K^ + '(X) from 
the left one K ( -~\X) by means of an isomorphism. 



3 The i^-matrices and open spin chains 

The integrability at the boundary is described in terms of two scaterring matrices K^'(X) |16j . 
The compatibility between bulk and boundary scattering lead us to an algebraic condition at one 
of the ends of an open chain, which reads 



R l2 {\-^K { c\\)R 2l {\ + ^)K { i\^ = K { 2 \fi)R 12 (\ + f i)K { 1 \X)R 21 (X-fi) (14) 

where K ( ~\X) = K^(X) ® I N and K { i\x) = In ® K^(X). 

When the bulk i?-matrix satisfies the properties fllQtil3[) one can follow a procedure devised 
in [PTl [18] to obtain the matrix K ( - + \X) at the opposite boundary At this point we note that 
property ffTSl plays the role of the crossing symmetry [18] and the matrix K^ + \X) is then fixed 
by the following isomorphism, 

*«(A)=lf(-)(^-A)«. (15) 

As a consequence of that we are left with the task to solve a single reflection equation. Here 
we shall be searching only for diagonal solutions of the reflection equation (I14p . namely 

N 

^ H W^ifH(A)e ¥ . (16) 

a=l 

The next step is to substitute the ansatz f[T6"j) in the reflection equation (fT4"l) and look for the 
simplest relations constraining the unknown matrix elements Ka \X). Among the many functional 
equations we find that there exists two of them that fix the ratios between the first two and the 
last two elements in a rather suitable way, 

K n {fx) fp^;^ - ^Rfi^ + a) + Rii(x - a*)^:?(a* + a) 

Rlj(X - ^)R 2 2 \^ + A) + ^^;i(A - fiSQiji + A) 



(19) 



and 

K { ~\[i) K^~\(X) N—l,N\ + N K N-1,NV A - N + tt N ,N-l\ A + W h n-i,nK* - N 

— r~\ = rA • (18) 

*n-i{V) R N _ 1 ' N (X + /i)^ atIi (A -fi)+ ' R N '_ 1N {\ + n)R N _{ N {\ - y) 

A jv_ii A ; 

Taking into account the explicit expressions for the Boltzmann weights we conclude that the 
right-hand sides of (I17II18P are independent of variable A only if the ratios K \_^ and K ^_ ) — 

K x (A) K N _ 1 (X) 

satisfy the following relations 

K t\X) _ sinh[A - A] K { -\\) _ sinh[^y_i - A] 
^(-)( A ) ~ sinh^ + A]' K^li(X) ~ sinh[/3j\r_i + A] ' 

where /?i and Pn-i are free continuous parameters. 

In the particular cases N = 2 and N = 3 the previous analysis already provides us a proposal 
for the K(~\\) matrix. The strategy for iV > 3 is to substitute Eq. (fr9|) into the reflection equation 
(1141) and to search for further relations that are able to determine the remaining ratios among 
next-neighbor amplitudes. By performing this analysis up to N = 5 we observed that some of the 
functional equations can be fullfiled once we extend the ansatz ffl9l) to any ratio ? 

^i!W = *4, „=!,...,*-!. (20) 
Kt\X) sinh[/3„ + A]' 

Now by substituting this proposal back to the reflection equation and after systematic algebraic 
manipulations, we find that the parameters j3 a have to satisfy the following recurrence relation 

&+i-& = ^r, a = l,...,iV-l. (21) 

Putting together all the above results we conclude that the most general solution for the 
amplitudes are, 

a-i , 1 N-l , 

Kt\X) = J] sinh[£_ + ,1 - ^(i - 6) - A] J] + ,1 - ^(i - 6) + A] (22) 

6=1 b=a 

where for later convenience we choose (3\ = £_ + + ^ such that the variable £_ is a free 
continuous parameter. We also emphasize that we have checked that the K-matrix f[2"2"j) satisfy 
the reflection equation f[T4"|) until N = 7. 



Considering the isomorphism (|T5l) one can easily derive that the form of the respective elements 
of K^ + >(X) are given by 

^ +) (A) = II sinh K+ + »? - l£G - 6) + A] II sinh ^+ + '? + ^5 + b) - X] (23) 



6=1 b=a 

where £+ is yet another free continuous variable. 

Having at hand the reflection matrices K^'(X) one can construct an integrable model with 
open boundaries following the double- row transfer matrix formulation proposed by Sklyanin |16j . 



f(A) =Tr A K%\x)T A (X)K ( A \X)T A (-Xy l 



where Xa(\) is the monodromy matrix of the corresponding closed chain with L sites, 



(24) 



T A (X) = RalWRal-iW • ..R M (X). 



(25) 



To obtain the respective Hamiltonian with open boundaries one needs to expand the double- 
row transfer matrix t(X) in powers A. The first derivative of t(X) is proportional to Tr A K A (0) 
which in our case is null for arbitrary values of 7 and N. In this situation we have to consider the 
expansion of t(X) up to the second order in the spectral parameter A. Considering that K A \X) 
has been normalized such that K A \o) = In, the expression for the Hamiltonian is [T9] , 

Hi, A 



H 



l^ H ^ + —dx K ^ W 

i=i 



+ 7 Tr A 

A=0 V> 



A=0 



1 



x=o PlA ~' 'MO) 



(26) 



where Hjj+i is the standard bulk Hamiltonian Hjj + i = Pjj+i -^Rjj+i(X)\ x=Q while ( is a constant 
proportional to the following identity matrix 



(In = Tr A 



d_ 

dX 



K A + \X) 



(27) 



Let us now present the explicit expressions of the open Hamiltonians in the simplest cases 
N = 2 and N = 3. For N = 2 we have the XX chain in the presence of both bulk and boundary 
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magnetic fields, 



L-l 



i=i 



where of, of and of are Pauli matrices. 



+ 



sin [7] 
2 

sin [7] 



cot 



cot 



7 7T 

s 2 4. 



7 7T 

! «+ + 2 + 4 



<r| (2f 



On the other hand for N = 3 we find that the corresponding open Hamiltonian up to an 
additive constant is, 

L-l 



1=1 



+ I sin[ 7 . fc + |] - 2 V / Sin[7]sin j 7 + £fcf] j [(fifS^ + SfSlJ S*S? +1 + 



+ 



+ 



sin [7] sin [7 + ^] 



2 sin[-<_ + I + f ] sin[-i£_ + | + f 



sm[-2<- + 7 + + ^ (^: 2 



y 



sin [7] sin [7 + 



irk 1 
3 -I 



2 sin[<+ + I + f ] sin[z£ + + \ + t> 



sin[2ze+ + 7 



27r^^ + ^3 ( ^ )2 



(29) 



3 J u 2 

such that S'f, Sf and 5"f denotes the standard SU(2) spin-1 matrices. 

In next section we shall consider the diagonalization of the double-row transfer matrix associ- 
ated to these models. 



4 Bethe ansatz analysis 

The purpose of this section is to present the eigenvalues A n (A) of the double-row transfer matrix, 

t(X) |0 n ) = A n (A) |0 n ) . (30) 

In the case of diagonal fT-matrices the diagonalization problem (13"U1) can be tackled within the 
algebraic Bethe ansatz framework. It turns out that the standard ferromagnetic highest vector 
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^o) is an exact eigenvector of t(X), 



)>=n®io> 4 , io) 4 



i=i 



/a 



V / 



(31) 



N 



playing the role of suitable reference state to start the Bethe ansatz analysis. 

An algebraic procedure to generate other eigenvectors of t(A) was first devised by Sklyanin 
for N = 2 [IB] . The central object in this approach is the matrix elements of the double-row 
monodromy matrix, 

r A {\) = T A (\)K { A \\)T A (\)- 1 (32) 
which also satisfies the reflection equation (fH]) . 

R 12 (\ - fi)T}(\)R 21 (X + rfTjip) = T»i? 12 (A + v)T}(\)R 21 {\ - ^ (33) 

The other eigenstates of t(A) were then constructed by exploring a set of commutation relations 
derived from the quadratic algebra (1331 . This algebraic analysis has been extended to tackle three- 
state vertex models [201 121] as well as the isotropic Heisenberg chain with arbitrary N up to the 
two-particle eigenstates [15J. In what follows we shall adapt the results of the latter work to include 
vertex models that are not invariant by charge symmetry. We shall not repeat the technical details 
discussed in [15] but only present the main relevant steps necessary to construct the eigenvalues 
A n (A). In order to describe that we shall first represent the double-mono dromy matrix as 



( % x {\) r li2 (A) 
r 2 ,i(\) r 2 , 2 (A) 



^i,jv(A) 

r 2iJV (A) 



(34) 



^ T N>1 (A) T7v i2 (A) ... T N:N (X) J 
Taking into account this representation, the diagonalization of the double-row transfer matrix 
t(A) becomes equivalent to the problem 



iY 



Y,k\ +) {\)% a {\) |</> n ) = A n (A) \4> n ) 



(35) 



i=l 



The first step in the algebraic framework is to reformulate (!35l) in terms of a suitable combi- 
nation of diagonal fields 7^ (A). It turns out that for arbitrary iV such linear combination is given 
by, 



M[t ) (2A)l , 



(36) 



where the j x j matrix M-^(A) is build up from the i?-matrix elements i?i,2(A) by the expression, 



M«(A) 



Rli(X) ijg(A) ... />'''';. | (A) Rl : i(X) 
Rll(X) m\) ... mi(X) R%{\) 



£g(A) Jig (A) ... /^,j(A) Jig (A) 



(37) 



Now by using Eqs fl35|36l) the eigenvalue problem can be rewritten as, 



N 



(38) 



i=l 



where functions w; + \A) are defined by 



w 



(+) 



" 'M«(2A) 



^; } (2A) 



^ (+) (A). 



(39) 



The basic property of the new diagonal fields ^(A) is that their action on the reference |0o) 
are always proportional to functions having a single power in L. More precisely we find that, 

, 2L 

(40) 



^(A)|0o) = ^ ) (A) r [ ff (A) |' J 0o) 



[p(A)M-A)r 



The functions -u^ ''(A) depend now on the amplitudes K\ } {X) by the expression, 



(-)/ 



i-l 



(-) 



^1,i(2A) I „ H 



(A) = ^" } (A)-E 



^r ; (A) 



(41) 



£ | iw£Vi(2A) 

where the second auxiliary j x j matrix M-J^A) is also build up from the J?-matrix elements 



.7 = 
r(~), 
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£i, 2 (A), 



<7 } (A) 



^'j(A) ... /*-/',.' (A) /^,.}(A) ^tJ(A) ... Bfi(\) 



■2,j+l 
'3+1,2 



;i(A) 

(A) R 



R i,i+i 



w(A) 

S53(A) 



^ji(A) 



?3,3+l 



.J',*+1 , 



(42) 



3 x 3 



To this point we have gathered the basic informations to provide a closed expression for the 
eigenvalue A (A) associated to the reference state \4>o)- In fact, by substituting Eq. (|40l) in Eq. (l38l) 
one finds, 



iV 



Ao(A) = ^m! +) (AK ( - ) (A). 



Ki(A)] 



2L 



(43) 



[P(A)P(-A)] J 

To obtain the other eigenvalues of t(X) we have to look for eigenstates generated by the action 
of the creation field %j(X) for i < j on the reference state |0o)- This analysis has been already 
carried out in [15] up to the two-particle states. Recall that in integrable theories the results for 
the two-particle state are enough to propose an educated expression for the eigenvalues A n (A). 
By adapting the conclusion of [T5] to our vertex model we find that the expression for the multi- 
particle eigenvalues are 

n rR i ' 1 mi 2i n 

= E r mw \^ w i +) (A) w !~ ) ( A) Yl Qi{\ Aj) (44) 
~i [p(A)p(-A)] fJi 

where Ai, • • • , A n are the variables parameterizing the n-particle state. 

The functions Qj(A,Aj) are determined in terms of the Boltzmann weights by the following 
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expression, 



#iSi(A-A*) #iJ+i(A-/i) 

^ ; : +i (a-^) ^ : :(A-/i) 



for % 



1 



^(A + zi) i?g(A + /i) 



^Ar-i,'i(A + /i) ^j?,T 1,2 (A + /^) 
-^iv-1,2 (A + A 1 ) -Rtv^ (A + A*) 



-^2jv(A - A 1 ) 



for 2 < i < JV - 1 



k Rn'UX - n) ^_i'i(A + (A + A 4 ) 



for i = N 



while the rapidities Aj satisfy the following Bethe ansatz equations, 



^2!i(Aj 



2L 



i +) (A,) ^(Ai) _ ^|i(2A i )^g(2A i ) - J R^(2A,) J R^(2A J ) " Q 2 (A„ A 



.(-)/ 



?2,2 



1.2/ 



>2,1/ 



^ +) (A,)^- J (A,) 



(-)/ 



Ki(2A,)]- 



n 



' Qi(Aj, Aj 



(45) 



(46) 



for j — 1, • • • , n,. 

We now have at hand the basic ingredients to exhibit explicit expressions for the eigenvalues 
A„(A). The first step is to simplify the expressions for w[ ± \X),--- ,w^(X) by using the K- 
matrices amplitudes (12211231) as well as the i?-matrix elements provided in section 2 and Appendix 
B. These simplifications require a considerable amount of algebraic work since we have to sum 
a number of distinct terms for each possible branch k. Fortunately, an analysis up to A^ = 5 is 
enough to exhibit the uniform dependence of these functions on N. We find that the final results 
factorize in terms of three types of products of trigonometric functions, namely 



a-l 



w { -\x) = n 



sinh[2A + ^(j-l)] 



L = Asinli[2A + i7+^(a+j-3)] 



Il sinh [£- 



A 



j 

JV-l 



i=i 



7 mk 

i ( 7 

2 N U 



2 n 



x 



I J sinh[£_ + X + t 



7 mk ,1 
2 ~ ~~N^2 



3)} 



(47) 
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and 



N-l 



n 



sinh[2A + ^(j-l)] 



L s inh[2A + z 7 + ^(a + j 



j 

N-l 



2)] 



smh[e + + A + z2-— (- 



-J)] 



x 



n sinh [^ 



A-.| + lr (iV+ 5 -i)] 



(48) 



The next step is to carry out similar algebraic simplifications for functions Qi(X, Xj). We notice 
that in order to make these polynomials as symmetrical as possible it is convenient to perform 
the shift Xi — > Xi — %\. Considering this change of variables and after some manipulations we find 
that the double-row transfer matrix eigenvalue are, 



An(A) 



fp(-A)" 


L 


L P(A) . 





a-1 



n 



sinh[A + '-£(} - 1)] 



X 



X 



N 

E 

o=l 
n 

J[ s inh[A - A* + ^ - (2 - a)^f] sinh[A - A 



2L 



L s inh[A + z 7 + ^(i-l)] 



sinh[A — Xi — i^) sinh[A — Xi + i 



w^{X)w^\X) 



iV J 



i I7rfc"| 
1 JV J 



n 



sinh[A + A, - %l\ sinh[A + A* + i| - if 



= | sinh[A + \ + %\ - (2 - a)*$] sinh[A + \ + %\ - (1 - a)*f] 
The corresponding Bethe ansatz equations for the shifted variables \ become, 



(49) 



/ sinh[A, + ^] \ sinh[e- + A, + g ] sinh[£ + - A, - g] 
VsinhfA, - z£] J sinh[£_ - A, + §£ ] sinh[£ + + A* - f$] 
_ " sinhjA, - Xj - g] sinh[A, + A, - g] 
y-sinh[A i -A i + ^]sinh[Ai + A i + ^] 

We conclude with the following remark. Note that the roots of unity branches on the left- 
hand side of Eq. (l50|) can be absorbed by performing the change of variables £_ — > £_ — and 
£ + — > £ + + After this transformation, the dependence of the Bethe ansatz equations on the 
roots of unity remains restricted to the two-body scattering amplitudes. 
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5 Conclusion 



The purpose of this paper was to solve the integrable vertex models based on the U q [SU (2)] algebra 
at roots of unity with open boundary conditions. We have solved the reflection equation and 
found one family of diagonal K-matrices having a free-parameter. For such diagonal boundary 
conditions we have been able to present the eigenvalues of the double-row transfer matrix and 
the corresponding Bethe ansatz equations. The next natural step would be to consider these 
vertex models with non-diagonal boundaries. In particular, to investigate if the functional relation 
approach developed for the open high spin XX Z quantum chain [221 12H [231 12S] can be applied 
to such roots of unity vertex models with non-diagonal i^-matrices. 



Appendix A: The projectors P 3 (^,k) 

For completeness we shall here present the explicit expressions for the projectors Pj(7,/c) given 
first by us in [T3]. In the Weyl basis these projectors can be written as decomposition of the braid 

S{<y,k), 

Pi(7,*) = n%^ (A.l) 
7=7 & ~ & 

where 

£. = (_i)i+i e *[5rCj-2)(j-D+27(j-i)]. (A.2) 
The corresponding matrix representation of the braid S'(7, k) is 

N 

S(j,k)= S c a j(7,fc)e M ®e a , c . (A.3) 

a>d;c>b 

where the amplitudes S^' b are given by, 



exp[ 2 -^(b- l)(d-l)+z 1 (b + d-2)} / ff(ffi, a -l)ff(ffi, c -l) 
HC-f,a-d) V H(f,<l-l)ff(f,6-l) 



/ g (27, a - l)ff (2 7 , c-l) x 
4/ g(2 7 ,d-l)g(2 7 ,6-l) a+h ' c+d ' ( j 
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n-1 

such that the auxiliary function H(X,n) = ]^[(1 — e^ x+ ~ 1 ^). 

1=0 

Altogether the above expressions provide us the explicit matrix expressions for the projectors 
for arbitrary and k coprime with N. 

Appendix B: Boltzmann weights for AT = 4 and N = 5 

In what follows we present the weights -R^(A) of the i?-matrix .R 12 (A) (Q for N = 4 and N = 5. 
We recall here that the only non-null amplitudes are those that satisfy the ice rule a + b = c + d. 



The non-trivial forty-four Boltzmann weigths for N = A are, 

iQ|}(A) = sinh[«7 + A] sinh[«7 + — + A] sinh[«7 + — + A] (B.l) 

R{'*(\) = R 2 2 j(X) = sinh[A] sinh[z 7 + ^ + A] sinh[«7 + ^ + A] (B.2) 

Rlj(\) = Rlj(\) = sinh[z 7 ] sinh[z 7 + ^ + A] smh[z 7 + ^ + A] (B.3) 

R{>*(X) = R 3 3 j(X) = sinh[A] sinh[A + ^] sinh[* 7 + ^ + A] (B.4) 

J#j(A) = R${\) = R 3 4(X) = Rti(X) = -|y s inhN sinh[* 7 + ^] 

x sinh[«7] sinh[«7 H — - — h A] (B-5) 

Rlf 3 (X) = Rlj(X) = sinh[z 7 ] sinh[^7 + ^] sinh[z 7 + ^ + A] (B.6) 

R{>l(X) = R 4 4 j(X) = sinh[A + ^] sinh[A + ^] sinh[A] (B.7) 



^m(A) = i^;S(A) = i#j(A) = i^(A) = -ej JsinhH sinh[z 7 + — ] sinh[A] sinh[A + ^] 

(B.8) 

J$4(A) = iSS(A) = i^;J(A) = Rli(X) = - efc ^ 8 inh[i 7 ] sinh[z 7 + ^] sinh[A] sinh[* 7 + ^] 

(B.9) 
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<i(A) 



(B.10) 
(B.H) 







rI%w - 


Rtiw 


— 


Rliw = 




X 


sinh[ry + 




— 


R^(X) = 




= 

X 


RltW = 

sinhfry — 
Rl : $(\) = 






sinh[«7 — 






it®A) = 


RIaW 




Rliw = 


R-ttW 




sinh[?7 — 



sinh[2ry] sinh[ - ] + sinh[A — ] sinh[A] 



sinh[A] 



(B.12) 



ITT k - — ITT k 

sinh[«7] sinh[^7 H — — ] + \J1e k sinh[A — ] sinh[A] 



ink. 



ink. 



(B.13) 
(B.14) 



>4,2, 



-,3,3/ 



ink. 



sinh[ry H — — ] sinh[«7 + — -] sinh[«7 — A] 
4 2 

• i r ■ i r ?7r ^n • i • i rn 

sinh[«7 H — —J smn^7 H J + smn|A —J smh|AJ 

4 2 4 

h[ry — A] sinhfry H — A] sinh[A] 



ink , , . , r ot/c 



(B.15) 
(B.16) 

(B.17) 

(B.18) 
(B.19) 
(B.20) 



where here we have = 1 or k — 3. 

For N = 5 the eighty-five Boltzmann weights are, 



i?M(A) = sinh[ry + A] sinhf^ + ^ + A] sinh[ry + + A] sinh[ry + + A] (B.21) 

5 5 5 

R±l{\) = R 2 2 j(\) = sinh[A] sinh[^7 + ^ + A] sinh[«7 + + A] sinh[«7 + + A] 



(B.22) 

2 1 / \ \ t-jI 2 / \ \ i r i ir ITlJv , r ^IttJv , r 327r/c -. 

Ri2\A) = R21W = Sinn [17 smh[?7 H — - — h AJ smh|ry H h AJ smh|ry H h AJ 

5 5 5 

(B.23) 
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R{>1(\) = Rl>{(\) =s inh[A]sinh[A+^]sinh[ i7 + ^ + A]sinh[z7+^ + A] (B.24) 



itg(A) = R%(\) = Rlj(\) = Rli(X) = -^2cosh[^] sinhN sinh[* 7 + ^] 

x sinh[A] sinh[ry + + A] sinhfry + + A] (B.25) 

5 5 

Rlj(X) = Rlj(X) = sinh[2 7 ] sinh[? 7 + ^] sinh[? 7 + + A] sinh[z 7 + + A] (B.26) 

o 



flM(A) = ^(A) = sinh[A]siiih[A + ^]siiih[A+^]siiih[i7+^ + A] (B.27) 

5 5 5 

<'4(A) = ^(A) = ^(A) = ^(A) = -sinh[A]sinh[A + ^]sinh^ 7 + ^ + A] 



x i / 2e fc sinh[z 7 ] sinh[z 7 H — ] cosh[— ] (B.28) 

V oo 

RliW = RltW = <'a(A) = Rli(X) = - sinh[A] sinh[* 7 + ^] sinh[* 7 + ^ + A] 

x W2£ fc sinhfry] sinh[ry H — ] cosh[^— ] (B.29) 

5 5 

i^(A) = i#*(A) = sinhH sinhk 7 + — ] sinh[ ?7 + — ] sinhk 7 + — + A] (B.30) 

5 5 5 

R\l(X) = 4^(A) = sinh[A]sinh[A + ^]sinh[A+^]sinh[A + ^] (B.31) 



i^;j(A) = Rl^(X) = i£g(A) = i^(A) = sinh[« 7 ] sinh[* 7 + 3^] sinh[A] (B.32) 

x sinh[A + ^] sinh[A + 2^] (B.33) 
5 5 

RliW = ^(A)=^(A) = ^(A) = -2sinh[A]sinh[A + ^] 

/ ■ n r i • n r ■ i r J^, . , r n CT ^i i r i7r ^i i ^^k 1 

x a / e fc smn[z 7 J smn^ 7 H — —J smn[z 7 + 2——\ smn|z 7 + 3— —J cosn[^^J cosn[ — - — J 
V 5 5 5 5 5 

(B.34) 



RliW = <2 5 (A) = Rtj(X) = Rti(X) = -\/e k sinhN sinh[* 7 + 3^] sinh[A] 

ITT k iTck 

x sinh[ry H ] sinhfry + 2 ] (B.35) 

5 ' 5 

Rl'l(X) = RH(X) = sinh[« 7 ] sinh[? 7 + ^] sinh[z 7 + 2^] sinh[^ 7 + 3^] (B.36) 

5 5 5 
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42(A) - 



sinh [27] sinh [27 + + sinh[A sinh[A] 

5 5 



. . r 2%nk . . r 3«7rA; 

Sinn [17 H h AJ sinh [27 H h A J 

5 5 

(B.37) 



- 42(A) 



4,3 (a) - 



4,4 (A) 



4,4 (A) 



43(A) 



4; 5 5 (a) 



41(A) 



■3,2/ 



2 sinh [17 smh|ry H —J cosh[^^J + smh|A —J sinh A 

5 5 5 



. _ r 3mk 
x sinh [27 H h AJ 



2,3/ 



smn[?7j sinh ^7 H — J + 2 smh|A —J smh|AJ cosh[^^J 

5 5 5 



x sinh fry H ] sinh [27 H h A] 

5 ' 5 

ink, 



ink, 



R 2 'i(X) = Rti(X) = sinh[A] sinh[A + — ] (sinh[A] sinh[A - — ] + 2e k sinh[z 7 ] 

5 5 

ink mk 
x sinh[ry + 3^^J cosh[^^J) 
5 5 



4,3(A) = 4J(A) = R tiW = 2 \l £ k sinh[?7 + ^] sinh[?7 



ink. 



ittJv ittJv itcJv 
x cosh[^^] sinh [A] (sinh [A] sinh [A — ] + sinh [27] sinh [27 + 3——]) 

5 5 5 

= i?4 2 (A) = sinh [27 H — — ] sinh[27 + 2^^](2e fc sinh[A] sinh[A — ] 

5 5 5 

x cosh[-— ] + sinh [27] sinh [27 + 3-—]) 

o o 

• 1 r 1 • 1 r Ink, . Ink, . . r n Ink, 

= sinh [27] sinh [27 H —J sinh [27 + 2——\ sinh [27 + o——\ 

5 ' 5 '" 5 

+ (sinh[A] sinh[— A]) 2 — 2 cosh[— — ] sinh[A] sinh[— — 

00 

. . r . . r ^I^k . /7rA; . Ink,. 

x (sinh [27] sinh [27 + 2 J + smn|ry H J smh|ry + 3 J) 

5 5 5 



A] 



i?5 2(A) = £fc sinh[«7 — A] sinh[A] sinh[A H ] sinh[A + 2"-^- 



ink. 



ink, 



= i?gJ(A) = 42(A) = ^(A) = 2e k sinh[? 7 - A] sinh[A] sinh[A + 

1 r ^k 1 L r 2mk, . ink, . . ^k, 

x cosh[— \ 2 cosh[^— J sinh [27 H — —J sinh [17 + o——\ 
5 V 5 5 5 



ink, 



4a(A) = 42( A ) = 45(A) = 2 sinh[? 7 - A] sinh[A] sinh [17 + 2—] 



>3,4/ 



>5,2, 



x cosh[ 



, ,2mk, . r «7r/c, . r ^TrA;, 

2 cosh[ J sinn|ry H J sinh [27 + 3 J 

5 5 5 



sinh [A] 
(B.38) 

(B.39) 
(B.40) 

(B.41) 
(B.42) 



(B.43) 
(B.44) 

(B.45) 

(B.46) 
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%7T k ITxl^j lift & 

B^'l(X) = R 2 5 i(X) = sinh[* 7 - A] sinh[? 7 + 3—] sinh[« 7 + 2—] sinh[z 7 + — ] (B.47) 

5 5 5 

Rl't(X) = iQ'f(A) = e k sinh[? 7 - A] sinh[A](sinh[« 7 + — ] sinh[2 7 + 3—] 

5 5 

x 2cosh[^]+sinh[A-^]sinh[A]) (B.48) 

R 4 / 4 (\) = R 3 / 3 (X) = sinh[* 7 - A] sinh[* 7 + 2—] (sinh[* 7 + — ] sinh[* 7 + 3—] 

^55 



+ 2cosh[— ]sinh[A- — ]sinh[A]) (B.49) 

i?3'g(A) = i?5 3(A) = 4e k cosh[— ] cosh[^^] sinh[A] sinh[A + — ] 

5 5 5 

ITT Aj 

x sinh[« 7 - A] sinh[? 7 + — A] (B.50) 

5 

RtiW = i$(A) = <'3(A) = i2K(A) = 2 cosh[^]^ cosh[^] sinh[* 7 + 3^] sinh[* 7 + 2^] 

ITT Jv 

x sinh[A] sinh[? 7 - A] sinhfry + — A] (B.51) 

5 

R 5 3 > 3 5 (X) = i2f;|(A) = sinh[* 7 + 3^] sinh[* 7 + 2^] sinh[* 7 - A] sinh[* 7 + ^ - A] (B.52) 

R\ 1(A) = (sinh 17 + 2— smh[? 7 + 3—] + Ae k cosh[— ] cosh — — 

5 5 5 5 

1TC k 17T Jv 

x sinh[A — ] sinh[A]) sinh[« 7 — A] sinh[« 7 H — A] (B.53) 

5 5 

R 4 / 5 (X) = Rli(X) = e k sinh[A] sinh[* 7 - A] sinh[* 7 + ^ - A] sinh[z 7 + 2^ - A] (B.54) 

5 5 



iq'g(A) = i2g'i(A) = sinh[« 7 + 3 ] sinh[« 7 - A] sinh[? 7 + A] sinh[« 7 + 2 A] 

5 5 5 



(B.55) 



,5 5/>- . , : . . , r 27r/c . n . , r 2wr/c , , . , r 3mk 



= sinh[z 7 — A] sinhfry H — A] sinh[ry H A] sinh[ry H A] (B.56) 

5 5 5 

where the possible values of k — 1, 2, 3, 4. 
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